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OPTIMIZATION OF TEMPERATURE DISTRIBUTION 
IN TISSUE BY MICROWAVE INDUCED HEATING POWER 


Piyanka DHAR,’ Rikhiya DHAR,” Ranajit DHAR* 


Abstract. An analytical study of optimal control. problem in distributed parameter 
system, described by one-dimensional single layered (homogeneous) bio-heat equation, 
is investigated so as to attain a beneficial therapeutic desired rise of temperature at a 
particular point of location of the tumour inside the tissue at the end of operation of the 
process by controlling optimally time dependent heating power induced by microwave 
and also by controlling surface cooling temperature. The heating power is constructed 
according to the well-known Beer’s Law [Karaa, Zhang and Yang, 2005]. The analytical 
investigation is carried out using “Maximal principle” [Pontrayagin et al., 1962] with a 
suitably constructed ‘Hamiltonian function’ by finite difference method. A numerical 
calculation of temperature distribution along the length of the tissue on various values of 
total time of operation of the process is computed. 
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Notations 
iC = specific heat of tissue, J/(kg °C) 
H = heat transfer coefficient between the skin and the ambient air, W m7/°C 
k = thermal conductivity of tissue, W m'/°C 
L = length of the tissue, m 
x] = position of tumour, m 
Xx = temperature, °C 
Xa = arterial temperature, °C 
Xo = initial temperature, °C 
uv = desired temperature to be attained, °C 
T = total time of the process, s 
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Q.(x,t) = spatial heating power induced by microwave, W m* 
u(t) = surface cooling temperature, °C 
Q,(t) = time dependent heating Power, Q in 
B = scattering coefficient, m" 
ty = switching time of operation of Q,(t), s 
ty = switching time of operation of u(t), s 
P = density of tissue, kg m™ 
re) = Dirac — delta function. 
@ = product of flow and heat capacity of blood, W m*/°C 
On = rate of metabolic heat generation, W m° 


1. Introduction 


In hyperthermia treatment, the temperature at the location of the tumour inside the 
tissue, described by bio-heat transfer equation, is raised to its therapeutic 
beneficial value while avoiding the damage of the healthy tissue due to 
overheating. In practice, microwave power level and cooling fluid temperature are 
the only variables which are accessible to direct control [Wagter, 1985]. It had 
been shown that surface cooling temperature could focus the microwave heating 
in deeper level of the tissue [Deng and Liu, 2002]. 


The fundamentals in the theory of optimal control of systems with distributed 
parameters were treated by [Butkovasky, 1969]. [Deng and Liu, 2002] obtained 
several analytical solutions to the bio-heat transfer problems with spatial heating 
on the skin surface using Green’s function method. Here one of the solutions was 
applied to study the temperature distribution in one-dimensional homogeneous 
tissue due to the application of spatial heating power constructed from well- 
known Beer’s law. [Dhar and Sinha, 1989] presented an analytical study on the 
strategy for achieving the desired rise of temperature of the tumour inside the 
tissue consisting of skin, fat, muscle and tumour layers by controlling surface 
cooling temperature with the finite difference method. 


A distributed optimal control problem, so as to attain a desired beneficial 
temperature of the tumour inside the tissue due to induced microwave, was 
analytically studied by [Dhar and Dhar, 2010] under conjugate gradient method. 
[Dhar and Dhar, 2011] carried out a theoretical study on optimal control problem 
in hyperthermia by controlling microwave induced heating power in multi-layered 
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tissue, where the ‘Maximal Principle’ was followed with a suitably constructed 
‘Hamiltonian’ by the use of finite difference method. The study of [Dhar, Dhar 
and Dhar, 2011] on the optimal distribution of time dependent point heating 
power, induced by conducting heating probe at the tumour site, was of 
considerable importance in hyperthermia treatment. 


[Karaa, Zhang and yang, 2005] presented a theoretical study on the temperature 
distribution in the tissue due to heating power induced by microwave, where the 
heat absorption in the muscle tissue was approximated by Beer’s law. 
[Kuznestsov, 2006]. studied analytically an optimal control problem in 
hyperthermia to attain the maximum rise of temperature on a located point of 
tumour at the end of operation of the process due to spatial heating. 


[Loulou and Scott, 2002] studied analytically a control problem on optimization 
of thermal dose in hyperthermia by using the conjugate gradient method. A 
fundamental study on optimization problems were thoroughly discussed on the 
concept of “Maximal Principle’ in [Pontrayagin et al., 1962]. 

[Wagter, 1985] presented a computer simulation model for optimization of 
microwave-induced hyperthermia. Later, a simulation-oriented optimization 


method to determine the input control variables in hyperthermia was obtained by 
[Wagter, 1986]. 


In case of spatial heating power, one of the important issues is to deal with the 
most typical one where the heat flux decays exponentially with the distance from 
the surface of the tissue [Deng and Liu, 2002]. Such spatial heating power 
induced by microwave was, in fact, constructed by well- known Beer’s law from 
which the heat distribution in the tissue could be well approximated [Karaa, 
Zhang and Yang, 2005]. 


In this paper, a distributed optimal control problem for a system described by 
Pennes bio-heat equation for a one-dimensional homogeneous tissue is 
analytically investigated such that a beneficial therapeutic desired temperature at a 


particular point of location of tumour x=, inside the tissue can be achieved at 
the end of total time of operation of the process by controlling optimally time 
dependant heating power induced by microwave Q,(t) [Wm] and also by 
controlling the surface cooling temperature u(t) [°C]. Here, the spatial heating 
power Q,(x,t)induced by microwave has been constructed according to well- 
known Beer’s law, given by, OQ, (x,t1)=Be™* Q,(t) where Q2(t) [W m7] signifies 
time dependent heating power and f is scattering coefficient [Deng and Liu, 
2002; Karaa, Zhang and Yang, 2005]. 
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In course of analytical investigation of this problem, it is seen that the optimal 
control variables Q,(t) and u(t) are singular controls. Thus, we have taken only 


two switching times ¢, and ¢,, for the sake of simplicity, in the consideration of 
distributions of Q,(t)and u(t) respectively. Here ¢, is the switching time during 
which the microwave induced heating power Q,(t) [Wm~] operates in the 
intervals (0,7,) and (¢,,7); t, represents the.switching time during which the 
surface cooling temperature u(t) is operative in the intervals (0,7,) and (¢,,7). 
Here T signifies the total time of operation of the process. 

The objective of this paper is to obtain the optimal values of controls Q,(t) and 
u(t) for different specified values of t, and t, respectively [Wagter, 1985] using 
“Maximal principle’ [Dhar and Sinha, 1989] with a suitably constructed 
‘Hamiltonian function’ by finite difference method. 

With these obtained values of Q,(t)and u(t), the temperature distribution along 


the length of the tissue at different times against various values of total time of 
operation of the process T are obtained numerically which display the therapeutic 
beneficial desired temperature on the point of location of the tumour inside the 
tissue. 


2. Mathematical analysis 


The one dimensional Pennes bio-heat equation can be written as [Dhar and Dhar, 
2011): 


Oy Ox 

| Ne 

Ot Ox 
Where according to Beer’s Law, Q,(x,t)=feOQ,(t); Q,(t) is the time 
dependent microwave induced heating power and / signifies scattering coefficient. 


+ ON, — H)t+Q,(x,1)+ Qm (1) 


Boundary condition: 


eh fy Duley} on x=0 (2) 
Ox 
1=%. onx=L (3) 
initial condition: 
HX, 0)= Xo (4) 


Let us reduce the system with distributed parameters, given by equations (1)-(3), 
to lumped parameters described by ordinary differential equations of independent 
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time variable, where finite difference method is used for dividing the segment 
(0, L) on the X -axis into p equal number of equal sub intervals of length h for 


L . : , ; 
h=—., and the temperature of the tissue at each point x = x, is designated by 
Pp 


quantities q,(t) ((=0,1, 2,... p) [Butkovosky, 1969; Dhar and Dhar, 1989]. 


y Tumour 


q» q p-1 Amar Am Qm- q q| {0 
Fig. A. 

Here, the temperature of the tumour at the point x = x, is represented by q(t), 

given in figure A. 


After discretising by finite difference method [Dhar and Sinha, 1989], the 
equation (1) stands as, 


Ch 
dt ° 
dq; 2 
=F, j =2,3,..,p-1 5 
| Ned P (5) 
The equations (2), (3), and (4) after discretising can be written as, 
1,74 
k= Ago uO} (6) 
qp > v4 (7) 
and g,O)=%, f=1,2,3,..., p (8) 


where 'F, and *F ; in equation (5) thus stand as, 


2 -Bh 
= 1 : k a me: k — k,Au(t) 7 Q,, 4 Be *"Q,(t) n b ‘ 
pe |h'(k+~Ah) h pch pch(k+Ah) pe pc pc 
(9) 


(10) 


k b Q,, , Be Q,() 
°F = ' =? + ; : 4 m +4 2 
J pch? (4. qj q;) (q, un) 


which are obtained with the help of equations (5)-(7). 


1 
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In order to attain a desired rise of temperature xv ata particular point of location 
of the tumour x = x, inside the tissue at the end of operation of the process T by 
controlling optimally microwave induced time dependent heating power Q,(t), 


[Wm] and also by controlling surface cooling temperature u(t) [°C], we would 
like to formulate the optimal control problem given as follows: to obtain the 
optimal control variables Q,(t), [Wm] and u(t) [°C] such that the quadratic 
deviation of the true temperature y(x,,7) from the desired temperature y° is 
minimized. 
Thus the objective function after discretising [Butkovosky, 1969; Dhar and Sinha, 
1989] stands, 

1 fees ” eee catonts 

D Ww —GA (T)} which is to be minimized. (11) 


Here, the temperature of the tumour on the point x =x, is represented by gq, (t) 


(given in figure 1) and q,,(7) designates the temperature on the point x= x, at 
the end of operation of the process T. 


Following [Butkovosky, 1969; Dhar and Sinha, 1989], the Hamiltonian can be 
constructed as 


p-l 
H=9,F,+>.9, °F, (12) 


j=2 
where @ ; (t) (j=1, 2,..., p—1) are adjoint functions. 
Let us write a functional J, given by, 


; vé d. p-l d ‘ 
Janse ~4n00) + [| Ste So, aha (13) 
. 


we will represent first variation of functional J as oJ. 


Thus, to obtain the optimality condition of control variables Q,(t) and u(t), we 
consider the stationary condition d/=0 for 6g AGF og 77) GL 2,64 p-D, 
oQ,(t) and du(t) by using calculus of variation and integrating by parts. 
This leads to a system of adjoint functions @ (¢), given by 

d aH 


— 9.()=—-—, (j=1, 2,..., p—1 14 
Fe, Ba, Gj p-l) (14) 


Optimization of Temperature Distribution in Tissue 
by Microwave Induced Heating Power 45 


With terminal condition 


p,7)=0, (Gj =1,2,....m—1, m+l,..., p—l) (15) 
?,T)=x% —4,(T) (16) 
The conditions of optimality of control variables Q,(t) [Wm] and u(t) [°C] 
fandins — (17) 
0Q, (t) 
and ae | =0 (18) 
Ou(t) 
Th 4 B Q <4 —jph 
us, ONE eg sien De ,(t) (19) 
Cc = 
Kx 
and t) =————— sign 0) (t 20 
u(t) Schuh) gn, (t) (20) 


From the equations (19) and (20) we note that Q,(f) and u(t) are singular 
controls. Thus Q,(t) assumes extreme values within the intervals (0,f,) and 
(t,,1) for switching time 1,. Similarly, the extreme values of u(t)is operative 


within the intervals (0, t,) and (t,,7) for switching time f,. 


Here it is assumed that the time control Q,(¢) [Wm] and u(t) [°C] are 
piecewise constant functions of time that change values at certain specified 
discrete instants considered as switching times [Wagter, 1985]. 

As ft, and ft, are switching times of control variables Q,(t) and u(t) respectively, 
the values of the optimal controls of Q,(t) and u(t) can be obtained from the 
equations (19) and (20). 


p-l 


Thus, ye "9, (t,) =0 (21) 


j=l 

and 2, (t,)=0 (22) 
It is to note that since we are dealing with the problem of optimal heating of the 
tissue, the control variable Q,(t) must assume the maximum value on the first 
time segment (0,7,) and consequently maintain minimum value on the second 
time segment (t,, 7) [Butkovosky, 1962]. 
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This aspect will make computer simulation one step forward towards finding the 
optimal values of control parameters Q,(f) . 


The value of optimal controls can be obtained with the help of equation (21) - (22) 
together with equations (8)-(22) by computer simulation [Dhar and Dhar, 2010]. 


3. Results and discussions 


Data used in computation are given as follows [Dhar and Dhar, 2010]: 


C = 3770Jkg"°C' 
p = 998kgm° 
K FO 5SW m’' °C" 
H A Wm? ec" 
f, = aye 
= ARSE 
L = 0.01 m, 
xy = 0.006 m 
o = 3000Wm* °C" 
B = 200m’ 
On = 33800Wm°* 
Xo = 25°C 
T = 600s, 800s, 1000s 


The computation has been carried out for various values of total time of operation 
of the process T and accordingly optimal microwave heating power Q2(t) and u(t) 
as well as the switching times f; and f during which they are operative have been 
changed [Dhar and Dhar, 2010]. 


In the figures we have mentioned, Q2(t) [W m”] as O(t) [W m”]. 


From figure 1, it is observed that temperature of the tissue [°C] attains desired 
temperature 43 °C at the point of tumour x = 0.006 m at the end of the process 
T = 600 s for 

Q(t) =1764 Wm”, 0 < t < 500s; Q(t) = 680 Wm” , 500s< t < 600s 


and = u(t) = 9.25°C, 0 < t < 4008; u(t) =6.25°C , 400 s<t < 600s. 
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Fig 1 : Temperature(°c) distribution of the tissue along the length of tissue(mm) for Q(t)=1764 wm”, 
0<=t<=500; Q(t)=680 wm”, 500<= t<=600 and u(t)=9.25°c, 0<=t<=400; u(t)=6.25°c, 400<= t<=600 
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Fig 2 : Temperature(°c) distribution of the tissue along the length of tissue(mm) for Q(t)=1525 wm”, 
O<=t<=600; Q(t)=300 wm”, 600<= t<=800 and u(t)=10.59°c, 0<=t<=500; u(t)=7.25°c, 500<=t<=800 
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In figure 2, the temperature of the tissue (°C) at the point of tumour x = 0.006 m 
rises to desired temperature 43 °C at the end of duration of the process T = 800 s 
for the application of 


Q(t) = 1525Wm”~ = 0 <t< 600s; 


O(t) = 300 Wm”, 600s<t < 800s; 
and u(t) =10.59°C,0<t< 500s; 
u(t) = 7.25°C, 500s<t < 800s . 


Temperature( °c) 


0 1 2 3 4 5 6 i 8 9 10 
Length of the tissue(mm) 


Fig 3 : Temperature(°c) distribution of the tissue along the length of tissue(mm) for Q(t)=1280 wm”, 
O<=t<=700; Q(t)=234 wm”, 700<= t<=1000 and u(t)=11.25°c, 0<=t<=600; u(t)=8.59°c, 600<=t<=1000 


Figure 3 displays the temperature of the tissue [°C] along its length [mm] due to 
application of 


Q(t) =1280 Wm” , 0 < t < 700s: O(t) = 234 Wm” , 700s<t < 1000s 
and 
u(t) =11.25°C, O0< t < 600s; u(t) =8.59°C, 600s<t <1000s . 


Here, the desired temperature 43 °C is attained on the tumour location x = 0.006 m 
at the end of operation of the forces T= 1000 s. 
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4. Conclusions 


It is observed that the temperature of tissue increases on the left side of the tumour 
at x =6 mm till it attains the beneficial desired temperature 43°C at the end of 
operation of the process and then the temperature of the tissue on the right side of 
the tumour decreases steadily to 37 °C (arterial temperature). Further, it is to note 
that as the total time of operation of the process increases from T=600s to 
1000 s, the first time segment of operation of the process (0,¢,) increases with 


the corresponding decrease of Q(t) [Wm] in this segment for the switching 
time t,. The surface cooling temperature u(t) [°C] increases in the first time 
segment of operation of the process (0, t,) as the total time of operation increases 
from T =600s to 1000 s for the switching time f,. 


Again, it is seen that the temperature of the healthy tissue on the both sides of the 
tumour are less than desired rise of temperature 43°C and thus the damage of the 
healthy tissue is avoided due to overheating. 


This analytical study of the optimal control problem may be used in case of 
computer- aided therapy planning in hyperthermia treatment. It can further be 
developed at different points of location of the tumour along with different length 
of tissue which may focus as a useful guideline to illustrate the versatility of the 
computer program. 
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